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Cross sections and form factors for e + e~ — > u>7T°, pi], and prj at center of mass energies of 
3.650, 3.686, and 3.773 GeV are measured using data samples collected with the BESII detector 
at the BEPC. Also, the branching fractions of tp(2S) — > ujtt , pi], and pi]' are determined to be 
(1.87to.62 ± 0.28) x 10 _s , (1.78t°;62 ± 0.17) x 10~ 5 , and (1.871^ ± 0.33) x 10 _s , respectively. 

PACS numbers: 13.25.Gv, 12.38.Qk, 13.40.Gp 



I. INTRODUCTION 

Form factors of the electromagnetic processes 
ip(2S) — > Loir (pr/, pr]') provide information on the 
strength of the electromagnetic amplitude in tp(2S) — > 
1~0~ decays. Such information is indispensable to 
separate the strong interaction amplitude for other 
ipi^S) — > 1~0~ decay modes, and even more impor- 
tantly, to determine the relative phase between the 
electromagnetic and strong amplitudes [1]. This phase 
has been found to be orthogonal in virtually all J/ip 
decays, such as 1+0" [2], l'O' [3, 4], 0"0- [5, 6], 
1~1~ [6], and NN [7] modes. It is suggested that 
this orthogonality also exists in ^(2S) decays and is 
universal [8]. 

The Born-order cross section for the electromag- 
netic process e+e~ — ► V + P (ujit , pr], pr]') in con- 
tinuum production can be expressed as [9] 

<JBorn(s) = ~^JY ■ \T V p{s)\ ' V V p{s), (1) 

where Vvp{s) = 3<Zyp, * s * ne momentum of ei- 
ther the Vector or Pseudoscalar meson in the VP de- 
cay, and the form factor |J r yp(s)| is an s dependent 
variable [10, 11]. 

By measuring the cross section of e+e~ — > VP in 
continuum production and correcting for the effect 
of initial state radiation, the form factor J-'yp(s) can 
be determined. In contrast, for e+e~ — > VP at the 
ip(2S) peak, the total cross section includes contribu- 
tions from both resonance production and decay and 
the continuum process, although the interference be- 
tween them can be neglected [9]. By separating the 
one-photon annihilation and ip(2S) resonance contri- 
butions, both the form factor at resonance and the 
branching fractions of ip(2S) — > VP can be obtained. 
The branching fractions of ip(2S) — > pi] and pi]' pro- 
vide useful information on the quark content of i] and 
i]' mesons [6]. 



In this paper, we report measurements for the cross 
sections and form factors of W7r°, pr], and pi]' at center 
of mass energies 3.650, 3.686, and 3.773 GeV, and the 
branching fractions of ip(2S) — > lott , pi], and pi]' at 
3.686 GeV. 



II. THE BESII DETECTOR 

The Beijing Spectrometer (BESII) is a conventional 
cylindrical magnetic detector that is described in de- 
tail in Ref. [12]. A 12-layer Vertex Chamber (VC) 
surrounding the beryllium beam pipe provides input 
to the event trigger, as well as coordinate informa- 
tion. A forty-layer main drift chamber (MDC) located 
just outside the VC yields precise measurements of 
charged particle trajectories with a solid angle cov- 
erage of 85% of 47r; it also provides ionization energy 
loss (dE/dx) measurements which are used for particle 
identification. Momentum resolution of 1.7%yl + p 2 
(p in GeV/c) and dE/dx resolution for hadron tracks 
of ~ 8% are obtained. An array of 48 scintillation 
counters surrounding the MDC measures the time of 
flight (TOF) of charged particles with a resolution of 
about 200 ps for hadrons. Outside the TOF coun- 
ters, a 12 radiation length, lead-gas barrel shower 
counter (BSC), operating in limited streamer mode, 
measures the energies of electrons and photons over 
80% of the total solid angle with an energy resolution 
of a E /E = 0.22/ VE (E in GeV). A solcnoidal mag- 
net outside the BSC provides a 0.4 T magnetic field 
in the central tracking region of the detector. Three 
double-layer muon counters instrument the magnet 
flux return and serve to identify muons with momen- 
tum greater than 500 MeV/c. They cover 68% of the 
total solid angle. 

In this analysis, a GEANT3 based Monte Carlo 
package (SIMBES) with detailed consideration of the 
detector performance (such as dead electronic chan- 
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ncls) is used. The consistency between data and 
Monte Carlo (MC) has been carefully checked in many 
high purity physics channels, and the agreement is rea- 
sonable. The generators EE2VP [13] and HOWL- 
VP [14], together with SIMBES, are used to deter- 
mine the detection efficiencies for the one photon an- 
nihilation and resonance decay processes. 

III. EVENT SELECTION 

The data samples used for this analysis consist of 
14.0 x 10 6 (1 ± 4%) ip(2S) events [15], 6.42 (1 ± 4%) 
pb -1 of continuum data at y/s — 3.650 GeV [16], and 
17.3 (1±3%) pb" 1 at the ?/>(3770) [17]. The chan- 
nels studied are lutt , pi], and pi]', where cu decays to 
7r + 7T~7r°, p to 7T+7T-, i] to 77, and rf to i]ir + ir~ . 

A neutral cluster is considered to be a photon candi- 
date if it is located within the BSC fiducial region, the 
energy deposited in the BSC is greater than 50 MeV, 
the first hit appears in the first 6 radiation lengths, 
and the angle between the cluster development direc- 
tion in the BSC and the photon emission direction 
from the beam interaction point (IP) is less than 37°. 
For the pr\ channel, tighter requirements are applied: 
the energy deposited in the BSC must be greater than 
100 MeV, and the angle in the xy plane (perpendicu- 
lar to beam direction) between the neutral cluster and 
the nearest charged track must be greater than 15°. 

Each charged track is required to be well fit by a 
three-dimensional helix, to originate from the IP re- 
gion, V xy = \Jv£ + Vy < 2 cm and \ V Z \ < 20 cm, and 
to have a polar angle | cos6*| < 0.8. Here V x , V y , and 
V z are the x, y, and z coordinates of the point of the 
closest approach of the track to the beam axis. 

The TOF and dE/dx measurements for each 
charged track are used to calculate Xpid^) vames 
and the corresponding confidence levels Probpiu{i) 
for the hypotheses that a track is a pion, kaon, or 
proton, where i (i = ir/K/p) is the particle type. For 
events involving the final state urn and pi]' , at least 
half of charged pion candidates in each event are re- 
quired to have Probpio^) larger than 0.01, while for 
events involving the final state pi], both charged pion 
candidates must satisfy this requirement. 

A. uj-n° channel 

For this channel, the events are required to have two 
good charged tracks with net charge zero and four or 
five photon candidates. A four constraint (4C) kine- 
matic fit under the hypothesis e + e~ — > 7r + 7r~7777 is 



performed. If the number of selected photons is larger 
than four, the fit is repeated using all possible com- 
binations of photons, and the one with the smallest 
X 2 is chosen. The confidence level of the 4C fit is re- 
quired to be larger than 0.01. In addition, we require 
that Xcom [18] for the assignment e+e~ — > 7r + 7r~7777 
must be smaller than that for e + e~ — ► K + K~j"/jj 
in order to suppress possible K + K^4 , y background. 

Among the four selected photons, there are three 
possible combinations to compose two 7r°s : (7172, 
7374), (7i73, 7274), and (7174, 7273)- Events where 
one and only one combination have both tt° candidates 
satisfying |M 77 - 0.135| < 0.05 GcV/c 2 arc kept for 
further analysis. 

Using the above selection, ir + ir~ 7r° 'tt° candidate 
events are obtained. Fig. 1 shows the 7r + 7r - 7r£ invari- 
ant mass distributions, where 7r° is the lower energy 
7T° from the two 7r° candidates. Clear lo signals are 
seen in all three data samples at yfs =3.650, 3.686, 
and 3.773 GeV. 

The 7T + 7r _ 7r° invariant mass distribution is fit with 
a shape for the signal determined by MC simulation 
plus a polynomial background, and 7.3TJ] 7, 31.2+^, 
and 8.6^3 3 events are observed in the data samples 
at Vs -3.650, 3.686, and 3.773 GeV, with 3.9c, 5.6<7, 
and 3.2cr statistical significance [19], respectively. 



B. pr\ channel 

For this channel, the events are required to have 
two good charged tracks with net charge zero and two 
photon candidates. A 4C kinematic fit to the hypoth- 
esis e + e~ — ► 7r + 7r~77 is performed, and its confidence 
level is required to be larger than 0.01 and larger than 
that for the assignment K + K~j-f to suppress possi- 
ble K + K~j~f background. The energy of each photon 
is required to be less than 1.7 GeV to reject back- 
grounds from ip(2S) — > 777'. An additional require- 
ment \m 7T+7T - - 0.776| < 0.15 GcV/c 2 is applied to 
further suppress background from non-p decay. Fig. 2 
shows the 77 invariant mass distributions for ,077 can- 
didate events. Clear i] signals are seen in all three data 
samples at </s =3.650, 3.686, and 3.773 GeV. 

The invariant mass spectrum is fitted with a 
shape for the 77 signal determined by MC simulation 
plus a polynomial background, and 2.3l 2 ;4, 29.2lg;g, 
and 5.81^6 even ts are observed in the data samples at 
Vs =3.650, 3.686, and 3.773 GeV, respectively. The 
statistical significances are 1.9cr, 5.5cr, and 3.3<r, re- 
spectively. 




FIG. 1: The ty + tv 7t" invariant mass distributions for 
e + e~ — > tt + it~ it ( Litx candidate events at yfs = (A) 3.650, 
(B) 3.686, and (C) 3.773 GeV. tt£ and ir° H are the low and 
high energy pions. The curves show the fit described in 
the text. 

C. prj' channel 

Here, events with four good charged tracks with 
net charge zero and two or three photon candidates 
are selected. A 4C kinematic fit is performed for the 
hypothesis 7r + 7r~7r + 7r~77. If the number of selected 
photons is larger than 2, the fit is repeated using all 
combinations of photons, and the one with the small- 
est x 2 is chosen. The confidence level of the 4C fit is 
required to be larger than 0.01 and larger than that 
for K + K ~ir + iT~-f-f. The two photons are required to 
come from r/ decay (|m 77 — 0.548| < 0.05 GeV/c 2 ). 
Background from tp(2S) — ► 7r + 7r~ J/tjj is rejected with 
the requirement that the mass recoiling from any 
ir+n- pair satisfies \mfj~ iX - 3.1| > 0.05 GeV/c 2 . 
In order to suppress background from non-p de- 
cay, an additional requirement \m^+^- — 0.776 1 < 
0.15 GeV/c 2 is applied, where m n + n - runs over all 
possible 7r + 7r~ pairs. Fig. 3 shows the rj-K + ir~ invari- 
ant mass distributions for p-n+Tt^rj candidate events. 
5.4+ 2 ' 2 events are obtained in the data sample at 
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FIG. 2: The 77 invariant mass distributions for e + e~ — ► 
/?77 candidate events at yfs — (A) 3.650, (B) 3.686, and 
(C) 3.773 GeV. The curves show the fit described in the 
text. 

\fs = 3.686 GeV by fitting the i]-k + tt~ invariant mass 
spectrum with an r/ 1 shape obtained by MC simulation 
plus a polynomial for background. The statistical sig- 
nificance for the rj' signal is 3.1cr. Only 1 event is ob- 
served in each of the two data samples at y/s =3.650 
and 3.773 GeV, and the respective background is zero 
and 0.64 events as estimated from sidebands; the cor- 
responding upper limit for the observed pr( event is 
calculated to be 4.4 and 3.9 events, respectively, with 
the scheme of J. Conrad et al. [20]. 

IV. SEPARATION OF CONTINUUM AND 
RESONANCE EVENTS 

At the ip(2S) peak, the cross section is from reso- 
nance decay, one-photon annihilation, and their inter- 
ference, all electromagnetic processes. The observed 
events are separated for the tp(2S) decay and one- 
photon annihilation processes using the scheme pro- 
posed by P. Wang et al. [9] , considering the detection- 
efficiency difference between the resonance decay and 
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FIG. 3: The rj-K + ir~ invariant mass distributions for 
e + e~ — > pn + -K~ri candidate events (solid histogram) and 
e + e~ — > prj Monte Carlo simulation (dotted histogram) 
with arbitrary normalization for y/s — (A) 3.650, (B) 
3.686, and (C) 3.773 GeV. The curve in (B) shows the 
fit described in the text. 

one-photon annihilation, and neglecting the tiny in- 
terference (see Table I), that is, 



N obs =C-(a 



e i>(2S) + a Cont. 



^Cont. 



where C is the integrated luminosity, a the cross sec- 
tion, and e the detection efficiency. Table I lists 
detailed information on the separation of the sig- 
nal events observed in the data sample at yfs = 
3.686 GeV. 

For the data sample at yfs = 3.650 GeV, the frac- 
tion of cross section from the ip(2S) resonance is neg- 
ligible (< 1.1 x 10~ 4 ) compared to that of the contin- 
uum, as estimated using the scheme proposed by P. 
Wang et al. [9]. Therefore we attribute all observed 
events to one-photon annihilation. 

For the data sample at \/s = 3.773 GeV , the reso- 
nance part of the cross section for e + e~ — > V(3770) — > 
VP (ljtt , prj, prj'), at Born order, is 

127rr r f 



(s - M 2 ) 2 + T 2 M 2 ' 



(2) 



where Tf is the partial width to the final state / 
(cj7t°, prj, prj') and is related to r ee and the corre- 
sponding form factor is [9]: 



a 3 



"ty(3770) 



\? 7 ("4(3770) ) I ' ''■ 



Since T t and T ee of V(3770) are 23.6 ± 2.7 MeV/c 2 
and 0.26 ±0.04 keV/c 2 , respectively [21], the fraction 
of the cross section from ^(3770) resonance produc- 
tion is negligible (w 2.1 x 10~ 5 ) compared to that 
from one-photon annihilation. The ip{2S) tail at 
\/s = 3.773 GeV contributes no more than 2.0% for 
the observed events in these final states. Therefore, 
the observed events are taken to be totally from one- 
photon annihilation. 



V. SYSTEMATIC ERROR 

Many sources of systematic error are considered. 
Systematic errors associated with the efficiency are 
determined by comparing J/ip and ip(2S) data with 
Monte Carlo simulations for very clean decay chan- 
nels, such as J/ip — > pir, W7T° and ip(2S) — ► Tr + ir~J/ip, 
which allows the determination of systematic errors 
associated with the MDC tracking efficiency, kine- 
matic fitting and photon selection efficiencies, etc. [22] 

To investigate possible background channels, we 
utilize three Monte Carlo samples : I) continuum 
channels from u, d, and s quark fragmentation gen- 
erated with JETSET7.4 [23]; II) ip(2S) -> anything 
generated with LUND-charm generator [24]; and III) 
DD pairs generated at \fs = 3.773 GeV. From these 
MC samples, we find that the background contami- 
nations arc negligible in the three data samples, ex- 
cept for the W7T° channel at \/s — 3.686 GeV. Events 
from ip(2S) — > bfir and cascade decays ip(2S) — > 
anything+J/tp, J/tp — > p°ir a can produce background 
in the M w mass region. The resulting estimate of con- 



tamination from these sources is Nbk g = 1.6 ±0.7, and 
the corresponding systematic error is (5.1±2.3)%. We 
take 7% as a conservative estimation and treat it as 
one source of systematic error for the urn channel at 
y/s = 3.686 GeV. 

Different s dependences of the form factor (- versus 
A-) result in changes of 2.5%, which is taken as one of 
the systematic errors. The uncertainties of the gen- 
erators, background shapes, and luminosities are also 
included. Table II lists all the sources of systematic 
errors, and the total systematic error is taken as the 
sum of the individual terms added in quadrature. 
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TABLE I: Cross section and event fractions for the ip(2S) production and decay and one-photon annihilation processes 
for oj7r°, prj, and prj' final states at yfs = 3.686 GeV. e is the detection efficiency determined from MC simulation. 
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TABLE II: Summary of systematic errors (%). Dl, D2, 
and D3 represent the data samples at yfs — 3.65, 3.686, 
and 3.773 GeV, respectively. 





U>7T° 


PV 


pv' 


Sample 


Dl 


D2 


D3 


Dl 


D2 


D3 


Dl 


D2 


D3 


Tracking 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


8.0 


8.0 


8.0 


Photon 


8.0 


8.0 


8.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Kine. fit. 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


4.0 


Bg. contam. 


0.0 


7.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Bg. shape 


13.6 


7.3 


11.6 


9.6 


4.5 


6.3 


0.0 


13.3 


0.0 


Luminosity 


4.0 


4.0 


3.0 


4.0 


4.0 


3.0 


4.0 


4.0 


3.0 


B(X -> Y) 


0.8 


0.8 


0.8 


0.7 


0.7 


0.7 


3.4 


3.4 


3.4 


MC fluct. 


1.4 


1.4 


1.4 


1.3 


1.3 


1.3 


2.1 


2.1 


2.1 


Generator 


1.3 


1.3 


1.3 


3.4 


3.4 


3.4 


2.8 


2.8 


2.8 


Fvp{s) 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


2.5 


Sum (a sys ) 


17.5 


15.0 


15.8 


13.3 


10.2 


10.8 


11.9 


17.9 


11.6 



VI. RESULTS AND DISCUSSION 



The Born order cross sections for 
calculated from 



X are 



(3) 



where N 



ClJ)S 



B(X^ Y)C-e MC ■ (1 + 5)' 

is the number of observed events 



in the final state from one photon annihilation, X is 
the intermediate state, Y is the final state, e is the 
detection efficiency obtained from the MC simulation, 
C is the integrated luminosity, and S is the radiative 
correction calculated according to Ref. [25] . With the 
cross section of e+e~ — > VP (ojtt , prj, prj') from Eq. 
(3) and using Eq. (1), we obtain \Typ\ at y/s =3.650, 
3.686, and 3.773 GeV. 

The branching fraction for tp(2S) — > X is calculated 
from 



B(ib(2S) -» X) 



• X- 



N, 



V>(2S) 



B(X -» Y) ■ e 



MC 



(4) 



where N°$ s) _ 



>x- 



is the number of observed events 



in the final state from it>(2S) decay. 

Tables III and IV summarize the results for e + e~ — > 
ivn°, prj, and pr]' at y/s =3.650, 3.686, and 3.773 GeV, 



and the branching fractions of ^>(2S) — > W7r°, ,077, and 
prj 1 at 3.686 GeV. The ratios of ip(2S) to J/ib branch- 
ing fractions are also listed in Table IV, where the J/rp 
branching fractions are taken from the PDG [21]. 
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FIG. 4: Energy dependence of the e + e~ — > 7* — > loty 
form factor. Curve (A) is calculated with Eq. (6), while 
curve (B) is calculated with Eq. (7). 

Fig. 4 shows the measured results of 1-7-^0 1 from 
our measurements, CMD-2 [26], and DM2 [27], and 

[4]: 



the calculated value of 1-7-^0 1 at s = m 2 /^ 



a (Pi \3 



r(v>(2S) -» ujtt 



(m 2 j^) 1 2 

i^o(o)i 2 " nui^^-nj/ip^p+p-) ' 

(5) 

where |.7w)(0)| = 2.3 GeV -1 [10], and the other 
quantities are taken from the PDG [21]. Curve (A) is 
predicted by J.-M. Gerard and G. Lopez Castro [11] 
as: 



|^W>(s — * 00 ) 



fujf-K 

3V2s 



(6) 



with the decay constants of u> and tt of f u = 17.05 ± 
0.28 and f„ = 0.1307 GeV. Curve (B) is predicted by 
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TABLE III: Cross sections and form factors measured for e + e — > luty", prj, and prj at yfs — 3.650, 3.686, and 3.773 
GeV. 



Channel 


Samples 


C (pb" 1 ) 


fjobs 


e(%) 


l + <5 


Co (pb) 


\Fvp\{ GeV" 1 ) 




3.650 GeV 


6.42 


7 o+3.3 
' •°-2.7 


5.09 


1.032 


24.3l* 1 ° ± 4.3 


o.o5it°:» 


L07T° 


3.686 GeV 


19.72 


17.3 + |^ 


4.98 


1.031 


19.2±f;f ± 2.9 


0.045lg : gg? 




3.773 GeV 


17.3 


8.6i 4 3 ;°3 


5.09 


1.028 


10.7±|"? ± 1.7 


034+ 008 




3.650 GeV 


6.42 


9 q+2.1 
z -°-1.4 


10.9 


1.028 


8.i^ : | ±1.1 


030+"' 014 
u.iwu_ 009 


PV 


3.686 GeV 


19.72 


16.0^;? 


10.9 


1.028 


18.4±f;| ± 1.9 


n rv/fc+0.011 

0-046t 0010 




3.773 GeV 


17.3 


5.8t 3 2 j 


10.7 


1.026 


7.8±3 5 ± 0.08 


030+ 009 




3.650 GeV 


6.42 


< 4.4 


4.33 


1.021 


< 89 


< 0.192 


pv' 


3.686 GeV 


19.72 


2-9±!:S 


4.43 


1.020 


18.61^3 ±3.6 


0.050l° : °- 




3.773 GeV 


17.3 


< 3.9 


4.56 


1.019 


< 28 


< 0.106 



TABLE IV: Branching fractions measured for ip{2S) — > cj7r°, pr;, and prj' . The corresponding J/i/) branching fractions [21] 
and the ratios Q h = B gtjJS^S are also given. 



Channel 


N i>(2S) 




e(%) 


B* ( as)^ (xlO- 5 ) 


Bj,^ (xlO- 4 ) 


Qn (%) 




1.4xl0 7 


13-9l!i 


5.97 


1-87±S:|| ± 0.28 


4.2±0.6 


44 +i.8 


PV 


1.4xl0 7 


13.21^° 


13.43 


1.781° j 7 ±0.17 


1.93±0.23 


9.2il:i 


pv' 


1.4xl0 7 


2.51?:? 


5.48 


l-87±i:ii ± 0.33 


1.05±0.18 


i7.8±i?:? 



Victor Chcrnyak [10]: 

|^7r"(s)| = |^7ro(0)| 



K- s )(M p 2 , 



(7) 



where m p and M p i are the masses of p(770) and 
p(1450), respectively. From Fig. 4, our results agree 
with the description of Eq. (6). 



VII. SUMMARY 

In conclusion, we determine branching fractions for 
ip(2S) — > lJ7t°, prj, and ,077' and the form factors T^a, 
T m , and T m , at Vs =3.650, 3.686, and 3.773 GeV for 
the first time. The branching fractions of i/j(2S) — ► 
o>7r° and pr7 in our measurement are consistent with 
those of CLEOc [28]. 
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